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ABSTRACT 

We report the discovery of a probable large scale structure composed of many galaxy clumps 
around the known twin clusters at z = 1.26 and z = 1.27 in the Lynx region. Our analysis 
is based on deep, panoramic, and multi-colour imaging (26.4' x 24.1') in V Ri' z'-bands with 
the Suprime-Cam on the 8.2 m Subaru telescope. This unique, deep and wide-field imaging 
data-set allows us for the first time to map out the galaxy distribution in the highest redshift 
supercluster known. We apply a photometric redshift technique to extract plausible cluster 
members at z ~ 1.27 down to i'=26.15 (5a) corresponding to ~ M* + 2.5 at this redshift. 
From the 2-D distribution of these photometrically selected galaxies, we newly identify seven 
candidates of galaxy groups or clusters where the surface density of red galaxies is signifi- 
cantly high (>5<r), in addition to the two known clusters. These candidates show clear red 
colour-magnitude sequences consistent with a passive evolution model, which suggests the 
existence of additional high density regions around the Lynx superclusters. 

Key words: galaxies: clusters: individual (RXJ 0848.6+4453, RXJ 0848.9+4452) — galax- 
ies: evolution — galaxies: formation — cosmology: large-scale structure of Universe 



1 INTRODUCTION 

Superclusters are the largest systems of galaxies, composed of mul- 
tiple clusters of galaxies, extending over 10-20 Mpc (e.g., Bahcall, 
Soneira 1984; Postman, Geller, Huchra 1988; Quintana et al. 1995; 
Small et al. 1998). According to the N-body simulations which 
successfully reproduce the observed filamentary structures at the 
local Universe (e.g., Peacock et al. 2001), clusters of galaxies at 
z ~ 1 are still in the process of formation (e.g., Moore et al. 1998). 
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Therefore, clusters of galaxies at this cosmological distance, espe- 
cially superclusters, are important sites where we can directly see 
the process of structure formation and evolution and mass assembly 
to cluster cores. 

Despite the importance in the cosmological context, however, 
only a few superclusters with the scale of > 10 Mpc have been iden- 
tified so far at high redshifts, especially beyond z ~ 1 (Connolly 
et al. 1996; Lubin et al. 2000; Tanaka et al. 2001). The limitation 
comes largely from the small field of view of large telescopes. 

The advent of the Suprime-Cam (Miyazaki et al. 2002), a huge 
format optical camera with a 30' field of view on the prime focus 
of the 8.2m Subaru telescope, has made it possible for us to view a 
>10 Mpc region all at once at high redshifts. Kodama et al. (2001) 
have taken this unique advantage of the Suprime-Cam on Subaru 
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Table 1. Log of the observations. 



Band 


Date 


Total Exp. 


lim. mag 

(5(7) 


Seeing 


V 


2000 Nov. 27 


96 min 


26.78 


~0"9 


R 


2000 Nov. 24 


90 min 


26.63 


~o;'9 


i' 


2000 Nov. 25 


30 min 


26.15 


~ o;'8 




2001 Mar. 21 


30 min 




~0"8 


z' 


2001 Mar. 22 


54 min 


24.88 


~ l'.'O 



to map out the large scale structure around the A85 1 cluster at z = 
0.41. They have found many filamentary structures and subclumps 
spreading over 5 Mpc from the dense cluster core, on the basis of 
photometric redshift (photo-z) analysis, an approach similar to that 
taken in this Letter. 

By utilizing the great light-collection power of the Subaru 
telescope, we now target the most distant supercluster ever firmly 
identified with spectroscopy: ie., the Lynx supercluster region at 
z ~ 1.27 composed of two known clusters, RXJ 0848.9+4452 at 
z = 1.26 and RXJ 0848.6+4453 at z = 1.27, which are spec- 
troscopically confirmed X-ray emitting clusters firstly reported by 
Rosati et al. (1999) and Stanford et al. (1997). The bolometric X- 
ray luminosities of RXJ 0848.9 and RXJ 0848.6 are 0.69^;^ X 
10 44 ergs s _1 and 3.3^0 5 X 10 44 ergs s _1 , respectively, which are 
derived with the Chandra observation (Stanford et al. 2001). 

Throughout this paper, we use the AB magnitude system 
(Oke, Gunn 1983). The adopted cosmological parameters are 
fin = 0.3 and Ao = 0.7, which gives a physical scale of 
8.37/ifQ 1 kpc arcsec -1 at the cluster redshift. 



2 OBSERVATIONS, REDUCTION, AND ANALYSIS 
2.1 Observations and Reduction 

The Lynx field was observed in November 2000-March 2001 in 
the V, R, i', and z' bands with the Suprime-Cam. Only eight of 
the ten CCDs of the Suprime-Cam were ready at the time of these 
observations. Consequently, each image has an FOV of 27' x 27' 
with a resolution of 0"2 per pixel. The log of the observations is 
given in TableQ 

The images are processed in a standard manner with IRAF and 
purpose-written software developed by us (Yagi et al. 2002). Flux 
calibration is performed with photometric standard stars from Lan- 
dolt (1992) for V- and i?-band, and spectrophotometric standard 
stars (Oke 1990; Bohlin, Colina, Finley 1995; Bohlin 1996) for i'~ 
and z'-band. We register the V, R, and i' images with the z' im- 
age, and match their PSFs to a fixed value, 1" (FWHM). Taking an 
overlapped region for all the images in the four bands, the final area 
used in our analysis is restricted in a 26'. 4 x 24'. 1 region, which 
corresponds to 13.3 x 12.1/if 2 Mpc 2 at z — 1.27. The regions at 
the edges of the combined image and regions around bright stars 
were masked and not used in the analysis. 

Since the i'-band image is the deepest among the four bands 
with respect to the passively evolving galaxies at the cluster redshift 
(i'(5<r)=26.15 mag or ~ M* + 2.5 atz — 1.27), we constructed an 
i'-band selected sample using SExtractor v.2.2.0 (Bertin, Arnouts 
1996). An area with more than 5 connected pixels with counts over 
2.8<T s k y is identified as an object. Photometry was then performed 
with a relatively small diameter aperture of 2" (16.7 'h^Q kpc) to 
keep the signal-to-noise ratios high enough for photo-z use. Galac- 



tic extinction is small, E(B — V) = 0.027 (Schlegel, Finkbeiner, 
Davis 1998), and hence neglected. After excluding 411 stellar-like 
objects, a total of 35143 objects, whose magnitudes are brighter 
than i' = 26.15 (5<r) in the ~ 590 arcmin 2 area, are contained in 
the final catalog. 

2.2 Foreground/Background Subtraction by Photometric 
Redshifts 

Our aim is to map out structures associated with the known twin 
clusters at z ~ 1.27. It is therefore required to remove unassociated 
galaxies in the foreground and background as much as possible, to 
maximize the contrast of the structures on the projected sky. 

Since spectroscopic measurements for the ~ 35, 000 galaxies 
are not practical, we exploit photo-z technique as an observation- 
ally efficient method to largely subtract the foreground/background 
populations. We input the VRi'z' magnitudes to the HYPERZ 
code (Bolzonella, Miralles, Pello 2000) to get estimated redshifts 
for all the individual galaxies in our i'-band selected sample. This 
code uses Bruzual & Chariot's (1993) stellar evolutionary code 
(GISSEL98) to build synthetic templates of galaxies for eight star 
formation histories: an instantaneous burst, a constant star-forming 
system, and six exponentially decaying SFRs with time-scales of 
from 1 to 30 Gyr. These models assume solar metallicity and the 
Miller-Scalo IMF (Miller, Scalo 1979), and internal reddening is 
considered using the Calzetti et al. (2000) model with Av varying 
between and 1.2 mag. 

It is ideal to have near-infrared magnitudes as well as opti- 
cal magnitudes to get better estimates for photo-z's for galaxies 
beyond redshift of unity, since the prominent spectral break fea- 
ture at 4000A for old galaxies starts to range out from the optical 
passbands at z > 1 (Kodama, Bell, Bower 1999; Connolly et al. 
1997). However, very limited FOVs of current near-infrared im- 
agers on large aperture telescopes (eg., 2' x 2' in the case of Sub- 
aru) make it impractical to do near-infrared imaging of such a huge 
area of our Suprime-Cam field (~ 30'). Nevertheless, by obtaining 
z'-band magnitudes at the longest optical wavelength, whose effec- 
tive wavelength is about 9200A, we can still deal with the galaxies 
out to z ~ 1.3 for photo-z estimates. Our target supercluster is 
supposed to be located at z ~ 1.27, hence the z'-band can still just 
catch the 4000A break features. In other words, this target is one of 
the most distant structures that can be explored by photo-z's based 
on the optical imagings alone hence with this huge field of view, 
except for the rest-UV-selected Lyman break galaxies at exotically 
high redshifts such as z > 3 (eg., Steidel et al. 1999). 

By using the best estimate of photo-z of individual galaxies, 
we select the only galaxies within the range of 1 ^ z p hot ^ 1-35 
to isolate the plausible members associated to the supercluster at 
z ~ 1.27. The number of these photometrically selected candi- 
dates for cluster members (hereafter photo-z selected candidates) 
is 2229. This rather broad range of redshift is intended to ensure 
that we include the bulk of the cluster population, allowing for the 
intrinsic errors on photo-z's. In fact, we find that 14 out of the 
16 known true members of the twin cluster, which are confirmed 
by spectroscopy (Stanford et al. 1997; Rosati et al. 1999), are as- 
signed photo-z's within the above range. We therefore estimate that 
the completeness of the members of our method is greater than 
80%. For blue cluster members, the estimated photo-z's tend to 
be smaller than the true values. It is probably because our photo-z 
technique misidentifies some blue cluster members as red galaxies 
at smaller redshifts (Kodama et al. 1999). For this reason, the dis- 
tribution of Zphot — Zgpcc is likely to be asymmetric around zero. In 
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Figure 1. Comparison between the photometric redshifts estimated with 
and without the if'-band data (ZyRizK and ZvRiz, respectively), for the 
galaxies in the ell, cl2 and gr3 down to K' AB < 23 (see text for the data). 
Panel (a) and (b) shows the red and blue galaxies, respectively, separated at 
the threshold colour (see text). 



fact, among the 14 spectroscopically confirmed galaxies, the stan- 
dard deviation of 5 galaxies with z p h ot > z spcc is 0.05, while that 
of 9 galaxies with z p hot < z spcc is 0. 17. Because of this, we take an 
asymmetric range for photometric members around z = 1.27 (ie., 
1 *S -Zphot ^ 1.35). The remaining 2 galaxies are assigned sig- 
nificantly lower redshifts around z ~ 0.55. Both of these galaxies 
have foreground objects near to their images, and the photometry 
of these member galaxies may be affected by these foreground ob- 
jects. 

To assess our photo-z's based on the optical data alone 
(Zvmz), we compare them to the photo-z's estimated WITH the 
if'-band data (Zvrizk), which are available for the selected areas 
of our fields. It is expected that the photo-z's at z ~ 1.27 with the 
if' -band have better accuracies since the if' -band neatly covers 
the longwards of the rest-frame 4000A break at this redshift. Here, 
we use the if'-band imaging data of the 3 regions (2' x 2' each; 
centered on ell, cl2 and gr3 in Figure[2j taken by ourselves with a 
near-infrared camera CISCO on the Subaru telescope The integra- 
tion time of each image is about 30 min and the limiting magnitude 
is about K' AB = 23 (5a). The comparison is shown in Figure Q 
The red and blue galaxies indicate the galaxies redder and bluer 
respectively than a threshold colour, (i' — z') — 1.05 — 0.008z', 
which corresponds to the typical colour of the local Sab galaxies 
(Fukugita, Shimasaku, Ichikawa 1995). The numbers of the red 
and the blue galaxies plotted in Figure Q are 48 and 200, respec- 
tively. For the red galaxies, we have good match between ZvRiz 
and ZvRizK, suggesting the accuracy of our optical photo-z's are 
reliable enough for these red populations. In fact, the standard devi- 
ation, average, and median of ZvRiz — ZvRizK for the red galaxies 
is 0.15, +0.06, and +0.04, respectively, if we exclude six highly 
discrepant galaxies with \ZvRiz — Zvrizk\ > 0.5. For the blue 
galaxies, however, we see a significant scatter, in particular, we tend 
to underestimate photo-z's if we lack the JC'-band. Therefore, we 
must be careful in treating the blue galaxies. In this paper, how- 
ever, we rely basically on the red galaxies, and hence this is not a 
big problem. 

Given the broad redshift cut for photometric membership, the 
above photo-z selected candidates should include significant frac- 
tion of field contamination which are not associated to the su- 
percluster. Later in § 3, we discuss the level of field contamina- 
tion in the structures newly found in this study. We stress, how- 
ever, this photometric selection has already removed the major- 
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Figure 2. The 2-D distribution of photo-z selected candidates on the sky. 
North is up, and east is to the left. The large and small dots indicate red and 
blue galaxies respectively. The contour levels are 2, 3, 4 and 5<r. See the 
text for details. The radii of large circles are O.lh^ Mpc at z = 1.27. 

ity (32914/35143^95%) of all the foreground/background galaxies 
while keeping the bulk (~80%) of the cluster members. This essen- 
tial process has now made it possible for us to unveil the large scale 
structure possibly associated to the known twin cluster and spread 
over the entire Suprime-Cam field, which are otherwise embedded 
in the foreground/background galaxies and invisible. 



3 LARGE SCALE CLUMPY STRUCTURE 

Figure [2] shows the spatial distribution on the sky of the photo-z 
selected candidates at 1 ^ -z p hot ^ 1.35. The large and small 
dots indicate red and blue galaxies respectively defined in §2.2. 
The numbers of the red and blue galaxies are 575 and 1654, re- 
spectively. The positions of the galaxies, X and Y, in the Figure|2| 
are expressed relative to the bottom-left corner of the Suprime- 
Cam image. The contours indicate the local surface density of the 
'red' photo-z selected candidates calculated from 10 nearest neigh- 
bors, corresponding to 2, 3, 4 and 5a. Here a corresponds to the 
scatter (standard deviation) of the local surface density in the case 
of the random distribution of 575 galaxies over the 590 arcmin 2 
field of view, ie., 0.97 arcmin -2 . We do not include the blue galax- 
ies here for the following reasons: (1) The accuracy of photo-z is 
relatively poor for the blue galaxies because of their intrinsically 
weak 4000A feature (§2.2). (2) The red galaxies trace clusters or 
groups of galaxies more neatly than the blue galaxies due to the 
morphology(colour)-density relation (Dressier 1980, Dressier et al. 
1997). In fact, many authors have exploited this technique and have 
successfully identified distant cluster candidates in their wide-field 
data (e.g., Gladders, Yee 2000; Goto et al. 2002). We note that the 
surface density of the blue galaxies is more than twice larger in the 
area of X < 4.5 (arcmin) compared to the other area, while the 
surface density of the red galaxies are nearly the same. This proba- 
bly comes from the fact that the CCD chips were different between 
these areas at the time of the observation; in the early phase, the 
Suprime-Cam consisted of two types of CCDs, SITe (ST-002) and 
MIT/LL (CCID20). The sensitivities of these 2 CCDs are different 
especially in the z'-band, where the limiting magnitudes differ by 
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Figure 3. The 2-D distribution of photo-z selected galaxies on the sky in 
the Subaru Deep Field. North is up, and east is to the left. The symbols and 
the contour levels are the same as in FigureEl See the text for details. 
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Figure 4. The j' — z' versus z' diagrams of the galaxies in the known 
clusters (ell and cl2) and our new cluster/group candidates (grl-gr7). See 
the text for details. 



~ 0.3 mag. Hence the number of detected galaxies at X < 4.5 with 
MIT/LL CCD's are increased compared to the other area. Since we 
count only red galaxies to define local surface density, however, the 
contours in Figure[2|are not affected by this effect. 

In Figure|2| we mark by large circles the areas where the local 
surface densities of the red galaxies are higher than 5a above the 
mean density. The ell and cl2 indicate the 2 known high redshift 
clusters at z ~ 1.27, RXJ 0848.9+4452 and RXJ 0848.6+4453, re- 
spectively. We newly find 7 candidates of galaxy clusters or groups, 
where the local surface density is as large as those of the known 
clusters (grl-gr7). 

As we noted in § 2.2, our photo-z selected member candidates 
should include significant amount of field contamination due to 
the generous criterion for photometric membership (1 ^ z p hot ^ 
1.35). However, the contrast of our new cluster/group candidates 
against the general field is always statistically significant. The num- 
ber density of photo-z selected candidates of red galaxies in the 
low density (< Oct) regions is 2.5 ± 0.2 h% Mpc~ 2 , while those 
of the high density regions within the circles of O.lh^Q Mpc radii 
centered on our cluster/group candidates including the known twin 
clusters range from 6.3±2.0 h% Mpc" 2 to 11.2±2.7 /i| Mpc" 2 . 
To verify the significance of the structures, we present the typi- 
cal galaxy distribution at z ~ 1.3 in the blank field for compari- 
son using the Subaru Deep Field data (SDF; e.g., Furusawa 2002; 
Kashikawa et al. 2003; Ouchi et al. 2003) in Figure [5] Here we 
have applied the photo-z using the same passbands (VRi'z') and 
have extracted only galaxies that fall within 1 ^ z p hot ^ 1.35, the 
same criteria that we used for the Lynx cluster field. The contours 
indicate the surface number density of the red galaxies within this 
redshift slice and the contour levels are the same as in Figure|2| As 
shown, only the two areas with >5<r peaks are found in the SDF, 
while 9 such regions are identified in the Lynx field. It is therefore 
likely that most of the structures that we found in the Lynx field are 
real. To further strengthen the reality of the large scale structures as- 
sociated to the twin cluster, we present below the colour-magnitude 
diagram, radial profile, and the estimate for cluster/group richness 
for each clump. 

Figure [4] shows the i — z' versus z' diagram for each clus- 



ter/group candidate including the two known clusters. Galaxies 
within a circle of 0,7h^ Mpc radius are plotted. The filled circles 
and the filled squares indicate the photo-z selected candidates of 
red and blue galaxies, respectively. The small dots show the other 
galaxies. For comparison, we also show the galaxies randomly se- 
lected from the low density (< Oct) regions in the bottom right 
panel of Figure|4| The number of galaxies shown in this panel is the 
expected number of galaxies within a circle of Q.lhjd Mpc radius 
calculated from the surface density of galaxies in the Lynx field. 
The solid line in each panel shows the red sequence of the cluster 
RDCS 1252 at z = 1.24 derived by Blakeslee et al. (2003). We 
applied an appropriate colour-term correction from HST/ACS fil- 
ters to Suprime-Cam ones. The dashed line in each panel indicates 
a predicted colour-magnitude relation at z — 1.27 for passively 
evolving galaxies formed at zy = 4.5 (Kodama, Arimoto 1997), al- 
though they are almost identical to the empirically fitted lines (solid 
lines). Not only in the two known clusters (ell and cl2) but also in 
our new candidates (grl-gr7), can we find many red galaxies con- 
sistent with the passively evolving galaxies which comprise the red 
sequences in these diagrams. 

We also show the radial profile of each cluster and clus- 
ter/group candidate in Figure [5] Only those galaxies brighter than 
i' = 26 are counted because of increasing incompleteness toward 
fainter magnitudes. Here we subtract the remaining field contam- 
ination using the SDF data (see above) from each bin after nor- 
malizing counts of field galaxies with referring areas. The aver- 
aged surface number densities of all galaxies and the red galax- 
ies are 7.8 ± 7.1 and 0.39 ± 0.65 h 2 o Mpc -2 , respectively, to 
the same magnitude limit of %' = 26. Given the large field of 
views, the field-to-field variation of the galaxy densities should 
be averaged over, if present, and hence the accuracy of statisti- 
cal field subtraction should not be largely affected by this effect. 
We see a pronounced increase in galaxy density towards the cen- 
ter of most of our cluster/group candidates. The only exception 
is gr6 that shows no excess of galaxy density at the center. This 
may be due to the misidentification of the center of the group, or 
intrinsically non-axisymmetric structure. We combine the profiles 
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Figure 5. The surface number density of all galaxies as a function of the 
distance from the cluster center. Cluster and group indicate the combined 
profiles of the two known clusters and those of our new candidates, respec- 
tively, and the best-fitted King's law profiles are also plotted by the solid 
lines. 



of known clusters and those of our new candidates (groups) sepa- 
rately, and find that both of these combined profiles can be fitted 
very well with the King's law profile, N(r) = iVn(l + r 2 / ' R 2 .)' 1 
(King 1966, 1972), with R c = 0.17 ± 0.02 h^ 1 Mpc (cluster) 
and R c — 0.23 ± 0.03 h^^ Mpc (group), which is consistent with 
estimations for local clusters (Bahcall 1975). 

Table[2|summarizes the properties of the known clusters and 
the cluster/group candidates. The J2000.0 coordinates of the each 
cluster/group candidate are given in columns 2 and 3. Columns 4 
and 5 show the surface number density of all the galaxies and the 
red galaxies only, respectively, within a circle of 0.7/i^q 1 Mpc ra- 
dius. We also only use galaxies at i' < 26. The remaining field 
contamination has been subtracted using the SDF data (see above), 
and the errors indicate the Poisson statistics based on the number 
of cluster galaxies and that of the subtracted field galaxies. In col- 
umn 6, we show the estimated richness of each cluster/group can- 
didate, using the No.5 indicator introduced by Hill & Lilly (1991). 
This indicates the number of galaxies within a 0.5 Mpc radius from 
the cluster center and within the magnitude range between mi and 
mi + 3 measured in the R band at 2 ~ 0.5. Since the R band 
at 2 ~ 0.5 roughly corresponds to the z'-band at 2 ~ 1.27, we 
can make a direct comparison with their measurements. The clus- 
ter richness estimated based on JV0.5 is 1-2 for ell and cl2, while 
0-1 for grl-gr7 (See Table 4 of Hill & Lilly 1991). Therefore we 
find that the candidates we newly find are poorer systems than the 
two known clusters. This is also suggested from deep X-ray data 
around the Lynx field taken with Chandra (Stanford et al. 2001; 
Stern et al. 2002). Using the Chandra image in the 0.5-2 keV band, 
we measure the flux of excess diffuse emission in the circular re- 
gion around each candidate with a radius of 35". The background 
(the unresolved X-ray background and non X-ray background) is 
estimated from a surrounding annular region. The result is shown 
in column 8 of Table [2] Note that gr6 and gr7 are outside of the 
Chandra field. We cannot detect any significant X-ray excess at 
more than the 3<r level for our cluster/group candidates except for 
the two known clusters. It should, however, be noted that the 3<r 
upper limit luminosity of each region (see column 9 of Table[2} is 



comparable to a typical luminosity of nearby groups of galaxies. 
Thus, the non-detection in the Chandra data does not necessarily 
imply that the regions we find have unusually low X-ray luminosi- 
ties compared with their optical richnesses. We argue therefore that 
our 'Suprime-Cam imaging + photo-z' method is a very powerful 
technique in finding clusters or groups in the distant Universe even 
beyond z > 1 out to z ~ 1.3. We note that this is practically the 
highest redshift structures of 'normal' galaxies that can be traced 
by optical images with our scheme (bracketing the 4000A-break 
region in the rest frame). 



4 SUMMARY 

We have presented the results of deep panoramic imaging of the 
Lynx supercluster field at 2 ~ 1.27 taken with the Suprime-Cam 
on Subaru Telescope. Our multicolour image covers an area of 
26' .4 x 24' . 1 (13.3 x 12. Mpc), allowing us for the first time to 
investigate large scale structure spreading around the known cluster 
regions at this high redshift. By applying the photometric-redshift 
technique, we have mapped out the spatial distribution of galaxies 
near the redshift of the supercluster down to ~ M* + 2.5. 

We have newly discovered seven cluster/group candidates. 
These candidates show red colour-magnitude sequences and cen- 
trally concentrated profiles, similar to those of the two known clus- 
ters in this field, indicating that they are likely to be real clus- 
ters/groups of galaxies around the Lynx supercluster, comprising 
a large-scale structure over ~13 /i^Mpc scales. The dumpiness 
of the candidates suggests that they are in early dynamical stages, 
being in the process of assembly to a massive cluster. 

The next key step will be to confirm with deep spectroscopic 
observations that the cluster/group candidates we discovered in this 
study are physically associated with one another and forming a 
real large-scale structure at z ~ 1.27. If confirmed, this will be 
the highest redshift structure of 'normal' galaxies ever identified, 
which will have a large impact on structure formation theories and 
will bring us invaluable information on the formation and evolution 
of clusters of galaxies at their early stages. 
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